Entangling operations and their implementation using a small amount of entanglement 
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We study when a physical operation can produce entangle- 
ment between two systems initially disentangled. The formal- 
ism we develop allows to show that one can perform certain 
non-local operations with unit probability by performing lo- 
cal measurement on states that are weakly entangled. 
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Much of the theoretical effort in Quantum Informa- 
tion Theory has been focused so far in characterizing 
and quantifying the entangement properties of multipar- 
ticles states. The reason for that lies, in part, in the 
fact that those states offer interesting applications in the 
fields of computation and communication. In practice, 
these states are created by some physical action (or oper- 
ation) involving the interaction between several systems. 
This suggests that the analysis of these operations with 
regard to the possibility of creating entanglement may 
play an important role in Quantum Information Theory. 
The first steps in this direction have been recently re- 
ported ||,|. There, given a Hamiltonian describing the 
interactions of two systems, it has been analyzed how to 
produce entanglement in an optimal way. 

In this letter we investigate which physical operations 
are capable of producing entanglement. This goal is 
partly motivated by the recent spectacular experimental 
progress in the field, where several physical set-ups have 
been recognized to produce entangled states ||. Thus, 
some of the questions we analyze in this paper can be 
stated as follows: given a machine acting on two sys- 
tems, can it create entanglement? If so, what kind of 
entanglement? The basic mathematical tool to answer 
these questions is the isomorphism introduced by Jami- 
olkowski |Q. We will extend such an isomorphism to re- 
late physical operations [equivalently, completely positive 
maps (CPM) £} on two systems and unnormalized states 
(positive operators E) acting on two other systems. This 
allows us to reduce the problem of the characterization 
of physical operations to the one of physical states. 

The relation between physical operations and states 
has a well defined physical meaning. In fact, from the 
isomorphism it follows naturally that given a physical 
operation £ acting on two separated systems A and B 
initially disentangled (but entangled locally to some other 
ancilla systems) we can always obtain the corresponding 
state E as an outcome. What is even more surprising is 
that, starting from the state E we can always perform 
some local measurements such that for certain outcomes 



the state of systems A and B changes exactly as if we 
had applied the corresponding operation £. 

This last property will allow us to answer an intrigu- 
ing question raised in the context of Quantum Informa- 
tion Theory. Let us assume that we have two qubits A 
and B at different locations and we want to apply some 
non-local operation. This situation raises, for example, 
in the context of distributed quantum computation [g], 
where non-local operations between different quantum 
computers are required. So far, it is known that one 
can use maximally entangled states, local operations and 
classical communication (LOCC) to perform that task as 
follows: we can teleport the state of A to the location 
of B, perform the operation locally, and then teleport 
the corresponding state back to A. In this process one 
has to consume two maximally entangled states (i.e. two 
ebits) apart from transmitting two classical bits in each 
direction jf|. However, it is known that for some kind of 
operations (like the controlled-NOT gate) one can econ- 
omize the resources, such that only one ebit is consumed 
. In fact, all the operations that have been studied so 
far J^-|p] require an integer number of ebits. We will show 
here that many operations require a non-integer number 
of ebits. In particular, if the operation can only entangle 
qubits weakly, the required number is much smaller than 
one. This automatically implies that many tasks in dis- 
tributed quantum computation can be performed with a 
much smaller entanglement than the one required so far. 

Let us consider two systems A and B at different 
locations, whose states are represented by vectors in 
the Hilbert space H.a,b, respectively, both of dimension 
d. Any physical action on those systems is represented 
mathematically by a completely positive linear map £ 
mapping the density operator p of those systems onto 
another positive operator £{p). The map can be written 
as 

£(P) = J2°KP°1 (!) 

k 

where Ok are operators acting on Ha <8> ~Hb- For the 
sake of generality, we have not imposed that the map 
preserves the trace of p, since we may be interested in 
physical actions that occur with certain probability [flo| . 

Our first goal is to determine when a given CPM is 
able to produce entangled states. Thus, we first recall the 
definition of separable operators. We say that a density 
operator p is separable with respect to systems A and B 
if it can be written as [llT|l 
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P = ^ \ a i)A{ai\ ® \bi) B (bi 



(2) 



for some integer n, and | a^A € Ha and \ bi)B € Ws- 
Otherwise we say that it is non-separable (or entangled) . 
Separable positive operators describe states that can be 
prepared using local operations and classical communi- 
cation out of product states, i.e. are useless for quan- 
tum information tasks that require entanglement. Dur- 
ing the last years, much theoretical effort has been de- 
voted to study the separability properties of operators 
jL2|. In particular, a necessary condition for separability 
of a given positive operator p is that p TjK > p3| . [l4[ . 
where Ta denotes transposition in TIa in a given or- 
thonormal basis Sa = This condition turns 

out to be sufficient as well when the sum of the dimen- 
sions of Ha,b does not exceed five (for example, for two 
qubits). In higher dimension there are examples of en- 
tangled states represented by non-separable operators 
whose partial transpose is positive For methods to 
study separability of operators which have positive par- 
tial transposition we refer the reader to |l^] . 

We can similarly define separable CPM; that is, £ is 
separable [Em if its action can be expressed in the form 



£(p) = ^(Ai ® Bi)p(Ai ® B t )\ 



(3) 



i=l 



for some integer n and where Ai and Bi are operators 
acting on Ha,b, respectively. Otherwise, we say that 
it is non-separable. Up to a proportionality constant, 
separable maps are those that can be implemented using 
local operations and classical communication only |fl7| , 
i.e. useless for several tasks in quantum information. 

From the defintions (0) and (||) it follows that if £ 
and p are separable, then £(p) is also separable. This 
reflects the fact that by local actions one cannot create 
entanglement. 

Let us consider a CPM £ acting on systems Ai and Bi. 
We define the operator Ea 1 a 2 ,b 1 b 2 acting on Ha ®Hb 
[where now Ha = Ha ± <& Ha 2 and Hb = Hb % <8 Hb 2 , 
and dim(HAi) = dim(HBi) = d\ as follows: 



Ea 1 a 2 ,b 1 b 2 —£(Pa 1 A; 
Here, P Ai a 2 = \<P)a 1 a 2 (<S>\ with 



P, 



B\B-2 



\Va 2 



(4) 



(5) 



and S — {\i)}f = i an orthonormal basis. In the definition 
([|) the map £ is understood to act as the identity on the 
operators acting on Ha 2 and H b 2 ■ The operator E has a 
clear interpretation since it is proportional to the density 
operator resulting from the operation £ on systems Ai 
and Bi when both of them are prepared in a maximally 
entangled state with two ancillary systems, respectively. 



On the other hand, we have 

£(PA 1 B 1 ) = d 4 tTA 2 A 3 B 2 B3(EA 1 A 2 ,B 1 B 2 PA3B3PA 2 A 3 PB 2 B3)- 

(6) 



This can be proved as follows. First, we can write 

d 2 tTA3B 3 (pA3B3PA 2 A 3 PB 2 B3) = PA 2 B 2 i 



(7) 



where T means transpose in the basis Sa 2 ® Sb 2 ■ Now, 
using (||) one can readily show that 



£(PA 1 B 1 ) = d 2 tTA 2 B 2 {EA 1 A 2 ,B 1 B 2 PA 2 B 2 )- 



(8) 



Equation (^) has a very simple interpretation. It reflects 
the fact that if we have the state Ea 1 a 2 .b 1 b 2 at our dis- 
posal, we can always produce the map £ on any state of 
systems A3 and B3 by performing a joint measurement 
locally such that both systems A 2 A 3 and B 2 B 3 are pro- 
jected onto the maximally entangled state (||). Of course, 
this will happen with certain probability. Below we will 
show how to implement CPM with unit probability using 
this method. 

The relations (^) and (^) induce a one-to-one corre- 
spondence between CPM acting on tensor product spaces 
and positive operators. In fact, this correspondece can be 
viewed as an extension of the isomorphism introduced by 
Jamiolkowski jf] to tensor product spaces. Using these 
relation it is very easy to show that: 

(i) £ is separable iff Ea 1 a 2 .b 1 b 2 is separable with re- 

spect to the systems (A1A2) and (B1B2). Thus, we 
can study the separability of CPM by studying the 
problem of separability of positive operators. This 
immediately implies that we can use all the results 
derived for the latter problem [Q . 

(ii) £ can create entanglement iff Ea 1 a 2 ,b 1 b 2 is non- 

separable with respect to the systems (A1A2) and 
(B1B2). In particular, we can always obtain a 
state whose density operator is proportional to 
Ea 1 a 2 ,b 1 b 1 out of separable states by entangling 
our systems locally with ancillas. 

(iii) Let us assume that E^j^ 2 B B > 0, where Ta 1 a 2 
denotes transposition with respect to Ai and A 2 
in the basis Sa- Then, if p TA i > we have that 
£(Pa 1 b 1 ) Tai > 0. If additionally Ea 1 a 2 ,b 1 b 2 is en- 
tangled (i.e. bound entangled), then we can al- 
ways produce bound entangled states out of non- 
entangled states by using the map £. We just have 
to entangle the systems locally with ancillas. 

(iv) If £ corresponds to a unitary action, the correspond- 
ing operator has rank one, i.e. it can be written as 
E= |*)(*|, where |W) e H Ai ®Ha 2 ®H Bi ®Hb 2 
is a normalized state. 
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Let us consider some simple examples with qubits 
(d = 2). First, let us assume that Ea 1 a 2 ,b 1 b 2 is an entan- 
gled state with positive partial transposition. According 
to (i) the corresponding completely positive map £ is 
nonseparable and according (iii) [£ (p)] TAl > for all p 
separable. But in this case, positive partial transposition 
is equivalent to separability and therefore £(p) 

is separable for all p separable. However, if we allow for 
input states that are locally entangled with ancillas, the 
final state will be (bound) entangled according to (ii). 

On the other hand, let us consider a family of phase 
gates of the form 

U(a N ) = e - MNCT * 1 ® a * 1 , a N ee tt/2 n , (9) 

where the cr's are Pauli operators. These gates are of the 
same kind as the ones used in the discrete Fourier trans- 
form The corresponding operator Ea 1 a 2 .b 1 b 2 = 
l^cwK^cwl) where 

\lp aN ) = COs(a N )\'i> + )A 1 A 2 \'i> + )B 1 B 2 

-ism(aj V )|* + )A 1 A 2 |* + )B 1 B 2 , (10) 

and |$+) and |^ + ) are Bell states. 

In the following, we will use the formalism introduced 
above to study how to perform non-local operations us- 
ing a small amount of entanglement. Let us consider 
a basis of maximally entangled states of systems A X A 2 
(and B]^) as =1(8 J7j|$), where Ui are a unitary 
operators and |<&) is defined in (g). If we perform a mea- 
surement in that basis and obtain |$i) J 4 lJ 4 2 and |$j)b 1 b 2 
the state of our systems will be £(Ui <g> UjPA 1 B 1 U- <g> Uj). 
Thus, we see that as a result of the measurement we ei- 
ther implement the CPM, £, or local operations followed 
by the CPM. Now we will show how to use this effect in 
order to perform non-local operations by using entangled 
states. We will restrict ourselves to the case of qubits, 
but our results can be easily generalized. 

Let us start considering the gates U{aj^) (||). The 
amount of entanglement of the corresponding state \tp arf ) 
( flO| ) is given by its entropy of entanglement 

E (^a N ) = ~Xl0g 2 (x) - (1 -X)l0g 2 (l - x), (11) 

where x = cos 2 (ajv) = cos 2 (7r/2 Ar ). On the one hand, 
E(ip a2 ) — 1, i.e. according to our discussion U(n/A) is 
capable of creating 1 ebit of entanglement. On the other 
hand, E(ip ai ) = 0, i.e. U(n/2) — —ia x (£> a x is a local 
gate. For N > 2, we have that E(i/j aN ) is monotonically 
decreasing with N. Note that for N sufficiently large, we 
can regard @ as an infinitesimal transformation and use 
the results of Ref. (2| to show that the gate can optimally 
create an entanglement proportional to ajy • We will show 
that in that limit C/(ajv) can be implemented with unit 
probability by using an average amount of entanglement 
also proportional to a at, assisted by classical communi- 
cation of approximately 2 bits in both directions. Thus, 
we provide examples of non-local gates which can be im- 
plemented using much less than 1 ebit of entanglement, 



the required entanglement being proportional to the en- 
tanglement capability of the non-local gate. 

We want to perform the gate on systems A3B3 and ob- 
tain the output state in systems A1B1. We assume that 
both systems A1A2 and B1B2 are in the state \ipa N ) an d 
we measure systems A2A3 and B2B3 in the Bell basis 
|*ii,i a ) = 1® <7ti,i 2 |*}) where ax,l = 1,01,2 = a x ,a 2 .\ = 
o~y , and 02,2 = o~ z . Note that all outcomes of the mea- 
surement are equally probable. If the outcome for A 2 A3 
is \^i lt i 2 ), we apply <Ji Xl i 2 to Ai and proceed analogously 
with B2B3. One can readily see that the resulting op- 
eration on A1B1 after this procedure will be: (i) U(ocn) 
if ii = jx] (ii) U{a N y = U(-a N ) if i\ ^ j\. Thus, 
with probability 1/2 we obtain the desired gate, whereas 
with probability 1/2 we apply U(—oti\r) instead, and so 
we fail. In order to apply the desired gate with probabil- 
ity one, we proceed as follows. If we fail, we repeat the 
procedure but with systems A1A2 and B1B2 prepared in 
the state \ip2a N )- With a probability 1/2 we will suc- 
ceed, and otherwise we will have applied U{— a at) 3 to 
the original state. We continue in the same vein; that 
is, in the fc-th step we use systems A1A2 and B1B2 in 
the state \4>2 k - 1 a N ) so that if we fail altogether we will 
have applied C/(-ajv) 2 _1 . For k — N we have that 

U{—(xn) 2 1 = —U(ctN), and therefore even if we fail 
we will have applied the right gate, so that the proce- 
dure ends. 

The total average entanglement which is consumed 
during this procedure is given by 

E[U(a N )]=J2U) E(^ aN _ k+1 )=a N f N , (12) 

k=l ^ ' 

where 

1 N 

f N = -J2 2 k E(^ ak )</«, = 5.97932. (13) 

fc=l 

In (|T2]), the weight factor of pk = (l/2) fc_1 gives the prob- 
ability that the fc-th step has to be performed. Thus, we 
obtain E[U(ajf)] < «jv/oo- Due to the fact that in each 
step of this procedure one bit of classical communication 
in each direction is necessary [ p0[ , the average amount of 
classical communication is given by 2 — (l/2) N ~ 2 bits. 

Although the procedure described above allows only 
to implement gates with "binary phases" = tt/2 n , 
any gate U (a) with arbitrary phase a can be approxi- 
mated with arbitrarily high accurancy by a sequence of 
gates of the form U{a at), consuming on average E < 
JooOt ebits of entanglement. Furthermore, this proce- 
dure allows to implement any arbitrary two-qubit uni- 
tary operation U. We can write any such operation as 
U = e~ lHt = lim„^ 00 (l — iHt/n) n , where H is a self- 
adjoint operator. We can thus apply infinitesimal gates 
U n = (1 — iHt/n) sequentially using an extension of 
the scheme described above. Note that after such an 
infinitesimal operation we can perform local operations 
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without consuming entanglement. This allows us to re- 
strict the form of the Hamiltonians to those that can be 
written as 



Hn = 



E 

fc— x,y,z 



Hi- 



(14) 



fc=i 



This can be viewed as follows. First, let us write H in 
terms of Pauli operators for systems A and B as 



H = a ■ a A + (3 ■ a B + a A ■ <ya B , 



(15) 



where 7 is a matrix, and a is the Pauli vector. If we apply 
an infinitesimal local transformation in A and B with 
Hamiltonians — a ■ a A and —0 ■ a B respectively, this will 
be equivalent to having H with a = /3 = 0. Moreover, 
prior to this operation and after the application of U n we 
can always perform local operations such that we obtain 
an evolution given by Hq (fi~4|), where the /i's are the 
singular values of Hq. Since the Hk commute, we have 
that the corresponding unitary operation is given by 



U„ 



^-iHxt/n ^-iHzt/n e ~ iH 3 t/n 



(16) 



a sequence of operations of the form (^), for which we 
already have provided a protocol. The required amount 
of entanglement is therefore given by E\j — /oo^(a*i + 
M2 + M3) ebits. 

Using the results of Ref . [§] , one can compare for small 
ol at (large iV) the average amount of entanglement used 
up to implement the gate (^J) with the maximal amount 
of entanglement which can be produced with help of a 
single application of the gate |2l|j. One finds that for 
that the ratio E[U(aN) /E create [J7(ajv)] is given 
by w 3.1268, i.e. the amount of entanglement required 
to perform the gate is about 3 times the amount of en- 
tanglement which can be created using this gate. Similar 
results are found for a general U — e~ lHt in the limit 

As we have restricted ourselves to single applications 
of the unitary operation, the protocol given here is very 
unlikely to be optimal in terms of the consumed entan- 
glement per application of the operation U. One might 
expect that in some asymptotic limit, the average amount 
of entanglement required to implement a gate U equals 
the amount of entanglement which can be produced us- 
ing this gate. For U{it/A) (^), we have that the protocol 
described above is optimal, as it consumes only 1 ebit 
of entanglement, which equals the maximal amount of 
entanglement which can be created with a single appli- 
cation of U(ir/4). Several other examples of this kind 
— all of them dealing with an integer number of ebits — 
have been proven to be optimal in [[^J^] . 

Finally, let us mention that we have restricted our- 
selves here to the implementation of non-local unitary 
operations. In fact, with the formalism introduce here 
one can extend the analysis to non-unitary operations 
and even to the implementation of non-local measure- 
ments. All these results indicate that the entanglement 



properties of a physical operation £ are directly related 
to the entanglement of the corresponding operator E. 
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In the procedure described above, only (TV — 2) steps 
have to be performed in practice, as it happens that the 
required operation to be implemented in step (TV — 1), 
U (7r/2), is a local operation which can be performed with 
certainty and without classical communication. 
As shown in [g| , this amount is strictly larger than the en- 
tanglement of the state associated to the operation, i.e. 
it is better to start with some initially entangled state 
rather than a product state in order to optimally in- 
crease the amount of entanglement. In the limit a at — > 0, 
the maximal amount of created entanglement is given by 
-^create [U{a N )] = 1.9123aiv. 
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